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We have investigated the constraints on the Cardassian model using the recent results from the 
Wilkinson Microwave Anisotropy Probe (WMAP) for the locations of the peaks of the Cosmic 
Microwave Background (CMB) anistropy spectrum. We find that the model is consistent with 
the recent observational data for certain range of the model parameter n and the cosmological 
parameters. We find that the Cardassian model is favoured compared to ACDM model, for a higher 
spectral index (n s w 1) together with lower value of Hubble parameter h (h < 0.71). But for smaller 
£f} , values of n s , both ACDM and Cardassian are equally favoured. Also, irrespective of Supernova 

constraint, CMB data alone predicts current acceleration of the universe in this model. We have 
f*"^ ' also studied the constraint on the as, the rms density fluctuations at 8h _1 Mpc scale. 
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PACS numbers: 98.80.Cq, 98.65.Es 



I. INTRODUCTION 



The recent era has been dedicated to precision cosmology. Different remarkable observations with great precision 
has set new direction in realising the details of the universal evolution. Over the last few years we have developed 
a firm broad outline about our universe to be spatially flat, homogeneous, isotropic on large scales and composing 
\ of radiation, ordinary matter including baryon and neutrinos, cold dark matter and a huge share of dark energy. 
£f) • There have been a host of phenomcnological models in support of this picture, which sets strong limits on the 
cosmological parameters. But the latest Wilkinson Microwave Anisotropy Probe (WMAP)p| has set a new challenge 
C*~) . for these models. WMAP has shown incredible precision in predicting different cosmological parameters including the 
positions of the peaks in CMB anisotropic spectrum. WAMP's accurate determination of angular power spectrum 
along with other astronomical data sets like SNIa observations 001 an d 2dF galaxy survey put significant limit on 
Qh' the cosmological parameters, which should be defended by the various phenomenological models. 

Since the significant observations regarding the SNIa Q,y|, several models both simple and complicated ones, have 
been used so far to justify the late time accelerating universe and at the same time proposing different ranges for 
the cosmological parameters. With ACDM among the simple and Quintessence^, K-essence0 and several others|j| 
as more complicated ones, each one aims in describing the universe in accordance to more precised observational 
agreement. One more thing common here is that these models accepts the idea of dark energy in any form of field or 
fluid. While there have been many other interesting ideas like the generalised Chaplygin gas[jj or rolling tachvons[loj. 
' a different proposal has been put by Freese and Lewis [Tl|. In this new approach the new cosmology is obtained 
by modifying the dynamics of the universe, namely, the Friedman equation, without seeking to another unknown 
- - 1 component like dark energy. 

They have proposed an universe composed of only radiation and matter (including baryon and cold dark matter) 
which expands in the speed up fashion if an empirical term, called Cardassian term, is added to the Friedman equation. 

H 2 = Ap + Bp n , (1) 

where A — and B and n are constants and are the parameters of the model. Here the energy density (p) contains 
only matter (p m ) and radiation (p r ), i.e, p — p m + p r . Since at present p m >> p r , p can be considered consisting of 
p m only at present. The new term, dominates only recently at redshift ~ 1. To provide the requisite acceleration of 
the universe as the outcome of the dominance of this term, n should be < 2/3. The model has two main parameters 
n and B. 

There are several interpretations for the origin of this new "Cardassian term" appearing in the Einstein's equation 
(1)E3- To accommodate flatness with only matter and radiation in this new picture, the critical energy density has 
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been modified and expressed as a function of the original one. Here the total density parameter at present Slo is 
defined as Pm0 + Pr0 where the redefined critical energy density p c is related to the actual one p ca {— 3Hq /8irG) as, 



Pc = Pc 



F(n,B) 



where 



D 



F{n,B) = [1+^^(1 
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O m o and Q r o are two parameters defined as fl m a 



and £l r o = ^ respectively. 



Now, as matter and radiation evolves, equation (1) can be expressed as the following, 

H 2 = A [p m0 a- 3 (l 
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From equation (3) it is very straight forward to express B in terms of fi r o and f2 m o to be 



B n-l 
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Substituting this expression in equation (4), equation (1) is finally recast in the following fashion 



H 2 — QmoHgd 4 




(4) 



(5) 



(6) 



One should note that the only model parameter appearing in the above expression is n. (For details discussion 
please refer ^3|) The other important point to notice is that the value n — corresponds to a ACDM model. This 
will be crucial when we shall talk later about the allowed region of the parameter space. 

In one of our very recent papers [l3]| we have considered the constraints imposed on this model's parameters from 
the CMB measurements like BOOMERanG and Archeop and the Supernova results. We found the model as a quite 
aspiring as an alternative to dark energy models with interesting limits on its parameter space (For other investigations 
in this regard please refer to 01)- Now in the light of the very precised measurements of WMAP it is certainly very 
much worthy to verify the consistency of the model and found the new restricted parameter space. We further study 
the constraints on in our model. 



II. THE CONSTRAINTS FROM WMAP 



The CMBR peaks arise from oscillation of the primeval plasma just before the universe becomes translucent. The 
oscillation of the tightly bound photon-baryon fluid is a result of the balance between the gravitational interaction 
and photon pressure and this oscillations gives rise to the peaks and troughs in the temperature anisotropic spectrum. 
Locations of the peaks at different angular momentum I, resulted from this oscillation, depend on the acoustic scale 
I A which in turn is related to the angular diameter distance to the last scattering D and the sound horizon at the last 
scattering si s as ^ . To a good approximation this ratio for I a can also be represented by the simple expression [T^ 

"A = 7T z~Z ■ (7) 

C s Tls 

where f(= J a~ l dt) is the conformal time and the subscript and Is represent the time at present and at the last 
scattering era respectively. c s is the average sound speed before last scattering defined as 

rT lB 

c s = fj^ I c s dr, (8) 

with 



i) 



C .7 2 = 3 + ^ (9) 

4 p r 



is the ratio of baryon to photon energy density. 
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FIG. 1: Contour plots of the first, second and third Doppler peak locations in the (fl m ,h) plane for different values of n and 
n s — 1. Full and dashed contours correspond to the WMAP bound of first and second peak locations and the dot-dashed 
contour corresponds to the Boomerang bound on third peak location. The region within the dotted line corresponds to the 
f — a confidence region for Q m for the Supernova data. The allowed region is the intersection of all the contours. 



To calculate the conformal time today and at the last scattering, we use eqn.(6), 



and 



to = 



df = 



df = 



1 



da 
X(a) 



da 



<0*> Jo *(«) 



(10) 
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where X(a) 



/( a+ iVo 



n m0 



Substituting the above expression in equation (9), we have the analytical expression for I a in case of this model 

r 1 da 

Jo 



1 71 

lA = — 



>0 X(a) 
f a ls da 
JO X(a) 



(12) 



In an ideal photon-baryon fluid model, the simple analytic relation for the position of the m-th peak and the acoustic 
scale is l m = m I a- But this simplicity gets disturbed by the different driving and dissipative effects which in turn 
induces a phase shift to the original position |15l|. This shift has been compensated by parametrizing the location of 
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FIG. 2: Same as Fig.l but with n a = 0.97. 



the peaks and troughs by 

l m = l A (m - <j> m ) = l A (m - (j> - 8<f> m ) (13) 

where <f> is the overall peak shift = <pi and 5<j) m = 4> m — <f> is the relative shift of the mth peak. This parametrisation 
can be used to extract information about the matter content of the universe before last scattering. Although it is 
certainly very difficult to derive analytical relation between cosmological parameters and phase shifts, Doran and 
Lilley 17] have given certain fitting formulae which makes life very simple. These formulae does not have a prior and 
crucially depends on cosmological parameters like spectral index (n s ), baryon density (uib = ilbh 2 ), Hubble parameter 
(h), ratio of radiation to matter at last scattering (n s ) and also fif s as representative of dominating present dark 
energy density at the time of recombination. We use these formulae (given in Appendix)to specify the positions of 
the peaks in our model and finally constrain the model with the most recent WMAP results. 

The locations of the first two acoustic peaks from the WMAP measurement Q of the CMB temperature power 
spectrum are 

l Pl = 220.1 ±0.8, 

1 P2 = 546 ±10; (14) 

notice that all uncertainties are within la. The location for the third peak is given by BOOMERanG measurements [Tsj 

l P3 = 825±{° 3 . (15) 

We have studied the locations of the first three acoustic peaks in the cosmological parameter space for the Cardassian 
model given by Eq. (6). The sole parameter of the model is n. The cosmological parameter space we have investigated 
is given by (r, n s , h, Wb, f2 m ). Throughout this paper, we have neglected the contribution from the spatial curvature 
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FIG. 3: Contour plot of the locations of first, second and third peak locations in (fi m ,n) plane for n s = 1,0.97 and h = 0.71. 
The full and dashed lines corresponds to the WMAP constraints on first and second peak locations and the dot-dashed line 
corresponds to the Boomerang constraint on the third peak location. The allowed region is the intersection of all the contours. 



and massive neutrinos, setting Qk = = 0. We have also neglected the contribution from the gravitational wave 
in the initial fluctuation. Because of the rather tight WMAP constraint on Wb, (w b = 0.0224 ± 0.0009 0), we have 
assumed Wb — 0.0224 in our subsequent calculations. We have also assumed the optical depth of the last scattering 
t = 0.11 which is within the range of the WMAP bound, r = 0.1661°;°^ 0. 

In Fig.l, we have plotted the contours of the first three acoustic peak locations, corresponding to the WMAP 
and BOOMERanG bounds given by Eqn.(14)-(15) in the (Q m ,h) parameter space for different values of n and for 
n s = 1.0. In the same figure we have also shown the 1 — a bound on fi m by fitting our model with the SCP (Supernova 
Cosmology Project) data Q- 

To obtain a fit to SCP data, we write the apparent magnitude as 

m(z,tt m ,n,M) = Uog w V l + M, (16) 

where A4 = M — 51og 10 i?o + 25 and Di = H^di is the dimensionless luminosity distance, given in ref The 
parameters f2 m ,« and M. are determined by minimizing 



X 



■> _ \ - [m exp (zi) ~ m(zj, n m ,n,M)] . . 



where of is the error in m exp (zi). While computing A% 2 = X 2 ~ Xmin contours we choose the minimum of M at each 
point of (Q m , n). With a Gaussian distribution, we put the 1 — a (correspoding to 68.3%) bound on f2 m for each set 
of n across x 2 min - 

In Fig. 2, we have plotted the same contours for the CMB peaks with n s = 0.97 and Supernova bounds. For both 
n s = 1 and n s — 0.97, there is no 1 — a confidence region of fi m for SCP data within the assumed range (0.1 to 0.4) 
for n > 0.4. 

In Table 1, we have given the the allowed region for Q, m and h coming from WMAP and BOOMERanG constraints 
on first three acoustic peak locations with and without the Supernova constraint. 

In Fig. 3, we have shown the same contours as in Fig.l and Fig. 2, but in (O m , n) plane fixing h = 0.71. 



III. CONSTRAINTS ON <j 8 



In this section we present the WMAP constraints on a% which is related to the amplitude of the galaxy fluctuations 
for different values of the model parameter n. erg, the rms density fluctuations averaged over 8/i _1 Mpc spheres, is 
determined by the COBE normalisation of the CMB power spectrum. This is done by essentially fixing the fluctuations 
at the last scattering for a given model. 

Doran et.al 19] have given an estimate of the CMB-normalized cg-value for a very general class of dark energy 
models just from the knowledge of their "background solution" \VL d (a), uj d (aj\ (here superscript "d" stands for dark 
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TABLE I: Limits on Q m and ft for different values of n s and n from CMB constraints 



Without Supernova 


n 


n a = 0.97 


n s = 1.00 




ft 




ft 





Q, m < 0.35 


ft > 0.65 


O ro <0.25 


ft >0.72 


0.1 


Q m < 0.38 


ft > 0.63 


fi m <0.29 


ft >0.68 


0.2 


0.13< On < 0.40 


0.60< ft < 0.78 


fi m <0.32 


ft >0.64 


0.3 


n m > o.i5 


ft <0.72 


0.13< fi m <0.35 


0.76> ft >0.61 


0.4 


fim > 0.18 


ft <0.66 


fi m >0.16 


ft <0.68 


0.5 






fi m >0.19 


ft <0.62 


With Supernova 




n s = 0.97 


n s = 1.00 


n 




ft 




ft 





0.21< fi m <0.35 


0.65< ft <0.76 


0.21< Q. m <0.25 


0.72< ft <0.77 


0.1 


0.17 < fi m <0.32 


0.66< ft < 0.77 


0.17< n m < 0.29 


0.68< ft < 0.78 


0.2 


0.13< n m < 0.27 


0.67< ft < 0.78 


n m < 0.27 


ft > 0.67 


0.3 


0.15< Q, m < 0.26 


0.65< ft < 0.72 


0.13< Q, m < 0.26 


0.66< ft < 0.76 



energy) and the cg-value of the ACDM model, with Oq = Oq(A): 



where f is the conformal age of the universe, a cq is the scale factor at matter radiation equality given by, 

fl r0 _ 4.31 x IP' 5 

aeq ~o m0 ~ h*(i-n d )- 1 yj 

and uj is the effective equation of state which is an average value for oj d during the time in which fl d is growing rapidly: 

1 = Jl at n d (a)/^(a)d\na 

As shown in one of our previous papers ^| , the extra term in equation 1 can be represented as dark energy component 
with slowly varying equation of state and in that case iv d o can be used as a fair approximation tow. In our case it is 
n-10. 

f2| f is an average value for the fraction of dark energy during the matter dominated era, before Q d starts growing 
rapidly 

/■lnatr 

Of = [lnotr - hiaeq]- 1 / n d {a) din a. (21) 

J In o oq 

For model without early quintessence, which is precisely our case, is zero. 

In order to compute a% for the ACDM model, we have used the standard definition 

2 _ f°° dK 2/ ,^ fSj^nry 1 



H 



a£ = -A'(«) -^-^ , (22) 



with r = 8h 1 Mpc and 

3+n s 



A"(«) = ^( — ) 2», (23) 
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TABLE II: Limits on erg for different values of n s and n and fi m — 0.25 



With Supernova 


n 


n s = 1 


n s = 0.97 









0.83< cr 8 <0.85 


0.77 < cr 8 <0.79 


0.1 


0.78< o-g <0.80 


0.72< a 8 < 0.74 


0.2 


0.73< a s < 0.75 


0.67< a 8 < 0.69 


0.3 


0.67< a 8 < 0.68 


0.61< as < 0.63 



T(k) is the matter transfer function describing the processing of the initial fluctuations, for which we use the result 



T(g)=M l + 2.34 g)[1 + 3- 



2.34q 



+ (16.1g) 2 + (5.46g) 3 + (6.71g) 



41-1/4 



with q = -ryMpc and T is [2C 



r = fi m /i exp 



-Slh 1 



2h 



5h is the density perturbation at horizon crossing and a fit to the four-year COBE data gives 
l0 5 S H (n s ,n m ) = lMn-°- 785 -° mlnn ™ exp[-0.95(n s - 1) - 0.169(n s - l) 2 ] 



(24) 



(25) 



(26) 



assuming that there is no gravitational waves and also no reionization. Taking in to account the reionization effects, 
the above expression has been corrected by Griffiths and Liddle j^: 



Sh ^ = 1 + 0.76t - 1.96t 2 + 1.46t 3 , 



(27) 



5h(t = 0) 

where r is the optical depth. The above correction is reliable upto r = 0.5. As mentioned earlier we have used 
t = 0.11 in our calculation. 

In Table II, we show the constraint on ug for different values of n and n s , normalizing the calculation at f2 m = 0.25. 

IV. DISCUSSION 



The first result is that the WMAP constraint on the first two acoustic peak locations together with the 
BOOMERanG bound on the third peak location restricts the model parameter n as n < 0.5 for n s = 1 and n < 0.4 
for n s = 0.97. One should note that to have current acceleration in this model, n < 0.66 is required, which means 
that CMB data alone predicts a late time accelerating universe. When one puts the Supernova constraint on top 
of this, the bound on n becomes n < 0.3 for both values of n s . Another interesting feature of our analysis is that 
in order to distinguish between a ACDM model and a Cardassian model, one has to impose a low value of h (e.g. 
h < 0.71) and higher value of n s (n s w 1). This has also been shown in Fig. 3 where we have plotted the contours of 
the peak locations in the n) plane fixing h — 0.71. It is clear that for n s — 1, the allowed region does not include 
n = which is precisely the ACDM case. But for n s = 0.97 it includes n = 0. In other words, for low values of h 
together with n s closer to 1, it appears that the ACDM model is more disfavoured than the Cardassian model. On 
the other hand for smaller values of n Sl ACDM and Cardassian model are equally favoured. An interesting point to 
note here is that WMAP Constraint on n s is n s = 0.99 ± 0.04. Hence although this remains still an unresolved issue, 
the fact that n s plays an important role in determining the nature of dark energy is one very interesting outcome of 
this investigation. We have shown it in the context of Cardassian model (the fact, that n s plays an important role 
in determining the nature of dark energy, has also been shown recently by Barreiro et.al |24| in the context of early 
quintessence model), but it would be very interesting to study this feature in a model independent way. 
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We have also calculated the WMAP bound on erg m our model for different values of n and n s normalizing the 
calculation at fi m — 0.25. Table II shows that erg gets lower for higher value of n and lower value of n s . Considering 
the fact that the joint WMAP/Large scale structure data set indicate a suppressed clustering power at small scale 
0, resulting a lower erg value, one can conclude that lower value of n s and higher value of n is favoured so far structure 
formation is concerned. 

The greatest challenge in cosmology today is to unravel the nature of dark energy. Although the WMAP observations 
has put stringent constraint on many cosmological parameters, but there remains still some uncertainty in some 
important parameters like n s , as, Q m which plays very crucial roles in determining the nature of the dark energy. 
Hence the future observations with stronger constraints on these parameters hold promise to uncover the nature of 
dark energy. 
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V. APPENDIX 

For completeness, we repeat here the formulas used in our search for parameter space. These fitting formulae are 
quoted from the cited literature |l7|. 

We assume the standard recombination history and define the redshift of decoupling zi s as the redshift at which 
the optical depth of Thompson scattering is unity. A useful fitting formula for z/ s is given by |23j : 

z u = 1048[1 + 0.00124«; b - - 738 ][l (28) 

where 

, 9l = QmSiw^- 2 ^! + 39.5i^- 763 ]-\ g 2 = 0.56[1 + 21. lu^- 81 ]" 1 , 

Wb = £lbh 2 and w m = fl m h 2 . 

And the ratio of radiation to matter at last scattering 

n s = Pr{zis)/p m {zis) = 0.0416UV 1 (z ls /10 3 ) . (29) 

The overall phase shift tp (which is basically the phase shift of the first peak) is parametrised by the following 
formula 



<p = (1.466 - 0.466n s ) a x r^ + 0.291^ , (30) 



where a% and ai are given by 



ai = 0.286 + 0.626 w b (31) 

a 2 = 0.1786- 6.308 w b + 174.9 w\- 1168 w\. (32) 

and ftf s is the average fraction of dark energy before last scattering which is zero for our case. 
The relative shift of the second peak (6 fa) is given by 

Sip 2 = c - cm a - c 2 rj~ C3 + 0.05 (n s - 1), (33) 

with 

c = -0.1 + (0.213 - 0.123Of s ) (34) 

x exp{- (52- 63.6f2f s ) w b } (35) 

ci = 0.063 exp{-3500 w 2 b } + 0.015 (36) 

c 2 = 6 x 10~ 6 + 0.137 (w fc - 0.07) 2 (37) 

c 3 = 0.8 + 2.3Slf s + (70- I26n? s ) w b . (38) 
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For the third peak, 

^3 = 10-dirJ+0.08(n,-l), (39) 



with 



di = 9.97+ (3.3- 3Qf s )w b (40) 
d 2 = 0.0016 - 0.0067fjf s + (0.196 - 0.22fif s ) w b 
(2.25 + 2.77Slf s ) x 10" 5 



w b 

The over all shifts for the second and the third peak is 4> + $4>2 and <f> + Sfo respectively. 



(41) 
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